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Abstract 

The critical behaviour of thin films containing quenched random impurities and inho- 
mogeneities is investigated by the renormalization-group method to the one-loop order 
within the framework of the n-component 4 -model.The finite-size crossover in impure 
films has been consdered on the basis of the fundamental relationship between the effec- 
tive dimensionality D e Q and the characteristic lengths of the system. The fixed points, 
their stability properties and the critical exponents are obtained and analyzed, using a 
e = (4 — _D e fj)-expansion near the effective spatial dimensionality D e Q of the fluctuation 
modes in D-dimensional hyperslabs with two types of quenched impurities: point -like 
impurities with short-range random correlations and extended (linear) impurities with 
infinite-range random correlations long the spatial direction of the small size. The dif- 
ference between the critical properties of infinite systems and films is demonstrated and 
investigated. A new critical exponent, describing the scaling properties of the thickness 
of films with extended impurities has been deduced and calculated. A special attention 
is paid to the critical behaviour of real impure films (D = 3). 
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1 Introduction 



The development of the theory of phase transitions in thin films in the last decades has 
been based mainly on statistical methods 1] such as, for example, the lattice mean-field 
(MF) approximation [21 E], the phenomenological Landau-Ginzburg (LG) approach jlj, 
the general scaling theory IE], arid renormalizat ion-group (RG) schemes of calcula- 
tion |7||H]. The general approach to the description of finite-size (FS) systems having 
different geometries, including the film geometry, is the Fisher FS scaling theory jSJ IE] • 
Within this framework the field-theoretical RG methods are applied to the investigation 
of the possible types of critical behaviour and the calculation of important characteris- 
tics like the critical exponents, scaling amplitudes and crossover functions [TJ The 
calculations are performed for systems with a general spatial dimensionality D. The 
predictions for real systems (D = 3) are made by extrapolations from the corresponding 
e = (Dy — .D)-expansions, where Du is the upper critical dimensionality [Oj). The method 
of e-expansions is widely and successfully applied to the investigation of the critical prop- 
erties of FS and infinite systems. It is proven to be particularly convenient for the study 
of complex models of real substances, where the competing effects play a substantial role. 

For FS systems one of the most interesting problems is the dimensional FS crossover (FSC) 
from the usual D-dimensional critical behaviour to the corresponding d = (D — D )- 
dimensional critical behaviour when the finite sizes Loi(i = 1,...,Dq) are less than the 
correlation length £ ~ (T — T c )~ u ; T c is the critical temperature jSJ [7j. As £ — > oo for 
T — > T c , the FSC (D — > d) in the asymptotic critical behaviour occurs always provided the 
number d of "infinite" dimensions Li ^> £, (i = 1, d) is larger than the lower borderline 
(critical) dimensionality Dl 

In this paper we shall study the FSC in D-dimensional hyperfilms (D — 1, d — D — 1) 
of thickness L with quenched random impurities and inhomogeneities (HIGH], which can 
be described by a properly chosen 4 -model of critical phenomena with an n-component 
fluctuation field (order parameter) 0(x) = {(f> a (x),a = l,...,n}. A special attention will 
be paid to the case of real films (D = 3). Thus we shall investigate a complex system 
with three competing effects: the thermal fluctuations, the FS, and the disorder. The 
main features of the critical behaviour in such complex systems can be described by the 
Wilson-Fisher RG recursion relations in the one-loop approximation 0. The one-loop 
approximation leads to a reliable prediction of the possible types of critical behaviour, 
namely, of the fixed points (FP) of the RG equations, their stability properties, and of 
the first corrections to the MF values of the critical exponents. We shall perform this 
programme for basic models of random impurities by introducing a new approach to the 
RG investigation of the FS effects. 

The RG investigations of films reveal the two limiting cases of the FSC: (i) thick (quasi- 
D-dimensional) films, where the ratio y = (L /£) tends to infinity (y ^> 1), and (ii) 
thin (quasi-<i-dimensional) films, where the ratio y tends to zero (|/< 1). In the former 
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case the upper borderline dimensionality is Djj = 4, whereas in the latter case this 
dimensionality is Djj = 5 [HIHj. The e-expansions for the cases (i) and (ii) are performed 
for e = (5 — D) — (4 — d) and e = (4 — D) — (3 — d), respectively. Here we shall present 
an unified RG investigation based on an e = (4 — /} e ff)-expansion, where D c g £ [d, D] 
is the effective dimensionality of the fluctuation modes (f>(x) ^T]. For this aim we shall 
apply the so-called integration over dimensionalities 8(— D g q — d) less than unity. This 
approach has been recently used to perform a 5 < 1-integration over the time axis for 
problems of the quantum critical behaviour in disordered superfluids ^2]- Working in a 
similar way one may obtain a description of the critical behaviour in impure thin films 
for any effective spatial dimensionality D c g £ [d, D] . 

The term effective spatial dimensionality (D e Q) has been recently introduced ^T] in an 
attempt for a description of the FSC as a smooth variation of the effective dimensionality 
of the fluctuation modes <f>(x) with the variations of the ratio y = (Lo/0- Although the 
precise dependence D e g(y) cannot be easily obtained, it seems intuitively obvious that in 
thin films D e g(y) ~ d for y <C 1 (the case of thin films) and D c g(y) ~ D = (d + 1) for 
y ^> 1 (the case of thick films or, equivalently, of almost infinite systems). 

The random impurities and inhomogeneities produce the effect of "random critical tem- 
perature" T c {x) which depends on the spatial vector x jHJ EH ES 14J. This type of 
(nonequilibrium, quenched) disorder is investigated by choosing a convenient, usually, 
Gaussian distribution for the random (nonequilibrium) temperature T c (x). The criti- 
cal behaviour of infinite systems both for point ^3] and extended ^3] impurities is well 
known; see also EI] ■ We consider the FS effect on the critical behaviour in films with 
these basic types of disorder. 

Within the model of point impurities with spatially isotropic short-range random corre- 
lations the FSC (D — > d) can be easily proven. Besides, we demonstrate that the critical 
exponents are smooth functions of D g q and, hence, of the ratio y = (L /£)- We shall show 
that the impure critical behaviour predicted by predicted by T. C. Lubensky JH] for in- 
finite systems with point impurities occurs in thin films for any effective dimensionality 
D g q > 2, i.e. for any ratio y = (L /£). 

The disorder of type extended impurities is described by infinitely-ranged random correla- 
tions along one or more spatial directions. For systems of a film geometry the appropriate 
choice of extended impurities is the case of one-dimensional (line) impurities orientated 
along the direction of the small size Lq and randomly distributed along the other spa- 
tial directions. This disorder is described b a modification of the model of short-range 
correlated point impurities in which the short-range correlations along the small size are 
substituted with infinite-range correlations. The length scale of the latters is much larger 
than both the correlation length £ and the thickness L$. So the strongly correlated along 
the small size point impurities behave like continuous uniform strings. In regards to the 
critical behaviour this disorder acts like point impurities with a short-range random dis- 
tribution along the large (infinite) dimensions Lj and an uniform distribution along the 
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small size L Q . The anisortopy of the random correlations lead to a quite unusual critical 
behaviour of the impure films. This behaviour is described in details. As a result of 
our analysis, we have identifies a new critical exponent describing the scaling law of the 
thickness Lq of the film. This exponent is analogous to the dynamical critical exponent z 
in disordered classical ^5] and quantum systems [To] . 

In Section 2 we present the model of consideration. The specific FS features of the 
model are discussed in Section 3. The integration in dimensionalities less than unity 
and the effective dimensionality D e g are introduced in Section 4. The RG investigations 
presented in Sections 5 and 6 rely upon the basic problems discussed in Section 3 and 
4. In Section 7 we summarize the results and discuss their applicability to real three- 
dimensional (3-D) films. The analogy between the FSC in thin films and the classical-to- 
quantum crossover [17] at low temperatures is also disenssed. 



2 Model 



We shall use the usual 4 -Hamiltonian (TC = H/T, ks = 1) of the n-component fluctua- 
tion field <fi(x). We shall stick to the usual notations [UJEI1EII in which the Hamiltonian 
7i can be written in the form 

H = \\ d D x{{V<f>) 2 + r(x)<f + 2u<p A ) . (1) 

In Eq. (JT|). u > 0. r(x) = r + ip(x) , where ip(x) is a random function intended to describe 
the disorder effects and r = ao(T — T cQ )/T c0 is the usual Landau parameter represented by 
the (bare) critical temperature T c o of the pure systems [<p(x) = 0]. The random function 
<f(x) obeys the Gaussian distribution [§1 ITU]: 

&(x)<f(x f )] R = A6(x-x f ) , (2) 

(A > 0). The function (p(x) is related to the "random" critical temperature T c0 (x) in the 
following way: T cQ (x) = T c0 [l — (p(x)/a ]. The true (renormalized) critical temperature 
T c will be a function of T c q and the renormalized values of the interaction parameters u 
and A. The averages (A)r of the physical quantities are defined by the functional 

integral 

(A) R = J l[V^(x)e-^S dD ^A &(x)} . (3) 

X 

We shall work with periodic boundary conditions. This means that we neglect the surface 
energy which is important in the investigation of other properties of thin films. This way 
of treatment gives the opportunity to investigate the net effect of the FS L on the critical 
behaviour of the impure films. 

The RG investigation is performed in the space of the wave vectors. Taking into account 
the lattice structure, the components q u of the ^-dimensional wave vector q = (q u ; v — 
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0, 1, ...,d) are confined in the first Brillouin zone (—ir/a < q v < ir/a). However, because 
of the long-wavelength approximation (LWLA), aq u <C tt, which is unavoidably included 
in field models like the Hamiltonian the upper cutoff A of the wave vector <f is much 
smaller than (ir/a), namelly, A = ^(ir/a), where 7 -C 1. While the LWLA approximation 
does not introduce restrictions on the dimensions Lj which are initially supposed to be 
almost infinite (Li 3> £ for any a < £ < 00), the same approximation restricts the 
variations of the thickness: a C Lo. Therefore our consideration is confined within the 
latter condition and our results could not be extended to exactly two-dimensional (2D) 
films (mono-atomic layers). It is convenient to write the D-dimensional wave vector q 
in the form q = (ko,k), where the wave number ko = (27rno/£o) with no = 0,±1,..., 
corresponds to the thickness L , and the ci-dimensional wave vector k with components 
ki = (2imi/Li), (rii = 0, ±1, ...,) corresponds to the (i-dimensional "subsystem" of the 
film. 

The short-range random correlations (J2J) correspond to point -like random impurities and 
inhomogeneities which are equally distrubuted along all D directions. In the space of the 
wave vectors q the Eq. (J2J takes the form 

[ V (q) V (q'] R = A5(q,q'), (4) 

where 8(x,y) is the Kronecker symbol. The disorder of type "extended impurities" is 
described by the Gaussian distribution given by 

[ v (fMf')] R = A5(f-f'), (5) 

where the vector r lays in the <i-dimensional subsystem; as given by the definition x by 
x = (xo,r). In the g-space, the Eq. (j3J) takes the form 

MkMk')} R = A6(k, k') , (6) 

or, equivalently, 

Wq)<p(q')U = Atf(£ k')6(0, k )5(0, %) . (7) 

These models of impure systems have been briefly explained in Section 1. Besides, they 
are described well in several preceding papers |TU1 IT5] and we shall not enter in more 
details. 



3 Film geometry and FSC 

In order to clarify the application of the Hamiltonian ((!} to FS systems we shall consider 
simple lattice sums which appear in the one-loop perturbation contributions to the " self- 
energy" parameter r and the interaction constant u: 

A m (r) = ^T, S m(k,r), (8) 
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where m = 1,2, k = \k\, V d = (L u ...,L d ) = V D /L Q , and 



1 i 

S m (k,r) = -Y: {ki + k2 + r)m . (9) 

In a close vicinity of the critical point T c , where r = £ -2 ~ 0, the upper cutoff A for the 
wave vector q can be ignored and the sums (JBJ) and © over k and ko can be extended 
to infinity because the essential contributions to these sums are given only by the small 
wave numbers (0 ~ q 2 ~ £~ 2 ). As we are interested in the critical behaviour in a close 
vicinity of the critical point, where £ 3> a, the essential contribution to the sums S m and 
A m will be given by the terms with wave numbers q v < £ and therefore we must choose 
the cutoff A£ > 1. This is consistent with the LWLA provided (a/£) < 7 <C 1. Under 
the latter conditions the cutoff A can be kept finite or set infinite without any effect on 
the results of the summation. As we shall apply the Wilson-Fisher RG method, it seems 
convenient to keep the cutoff A for the wave vectors k and neglect the upper cutoff for the 
wave numbers ko- It is important to keep in mind that our investigation of the critical 
behaviour of slabs is valid only for 

£A > 1 , L > a . (10) 



While the fc-summation in Eq. (jHj) can be always replaced by a fc-integration, the sum- 
mation over the wave vector component k can be transformed to an integration 
only if certain conditions are satisfied. The latter become clear from the result of the 
summation (j^J) over k G (—00, 00) for m = 1: 

where y(k) = L (k 2 + r) 1 ^ 2 . For y(k) <C 1, Si from (fTT]) coincides with the (ko = 0)-term 
in 0. In this case, the .D-dimensional film has a d = (D — 1) dimensional behaviour. 
The quantity y(k) is much less than unity for any k only if (L A) 1 and y(0) = y = 
(L /£) 1. These conditions are consistent with (fTUj) provided a <^ L <^ (1/A). The 
condition (LqA) <C 1 is satisfied in a broad interval of values Lq > a because of the fact 
that A (ir /a). Moreover, having in mind that the small values of the wave number k 
yield the essential contribution in the sums A m one may consider the condition L ^ £ 
instead of (LqA) 1. This new condition for y(k) <C 1 is weaker than (LqA) ^ 1. 
Therefore, the limiting case y(k) ~ does not necessarily correspond to an exactly two- 
dimensional (2D) film (a single atomic layer; L = a). Rather, we have shown that the 
exactly 2D films are beyond the scope of the LWLA and that one of the main limiting 
cases of our consideration is the quasi- 2D film defined by a Lq (1/A). 

The second limiting case is given by y(k) ^> 1. The cases y(k) 1 and y(k) ^> 1 are 
often referred to as the limits y(k) — > and y(k) — > 00, respecticely. For y(k) — > 00, 
the Eq. (JTTJ) yields a result that can be obtained by replacing the sum Si by an integral 
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over k 0} namely, by taking the continuum limit along the small size L Q . This corresponds 
to a .D-dimensional behaviour of the system when all summations over the wave vector 
components q u can be substituted with an integration over the vector q. The quantity 
y(k) tends to infinity for any k only if y — > oo. Thus our film is quasi-3-D for all Lq 3> £. 

The behaviour of the sum Si in the two limiting cases y(k) — > and y(k) — > oo exhibits 
the dimensional FSC: d — *■ D — (d + 1). Having in mind that A m+ i = —(dA m /dr), 
m = 1, 2, we see that this dimensional crossover is a property of all perturbation terms 
and, hence of the system as a whole. 

The present discussion is particularly important for RG studies of the FSC, where the 
perturbation integrals are calculated at r = 0. In such studies the asymptotics of the 
corresponding integrands are considered and, therefore, one must be sure that the product 
kL Q can reach the limiting values (kL ) -C 1 and (kL ) 1 within the limitations ()10|) [8 . 
The RG integrals are given by taking the continuum limit for the summation over k in (jH|) 
and by confining the integration in the limits k £ [A/6, A], where b > 1 is the RG 
rescaling number. For such integrals with a lower cutoff A/b the condition (AL ) ^> b 
for thick (almost-3.D) films at the critical point (r = 0) is satisfied for L ^> (60/771*). 
The condition for thin (quasi- 2D) films is (ALq) <C 1 and, together with a <ti Lq, we 
have a < Lo ^ ( a /7 7r )- All these conditions are consistent within the LWLA and the 
requirement for a criticality as given by inequalities JTUJ). 

Therefore the RG investigations are consistent with both the aims and the model chosen in 
this paper. The behaviour of the exactly 2D films is beyond the scope of our consideration 
but we can reliably investigate the FSC from quasi-3-D to quasi-2D films. 

4 Integration in noninteger dimensionalities 

An obvious disadvantage of all existing descriptions of the FSC and other dimensional 
crossovers is that the limiting cases can be easily proven and described but the 
intermediate case (0 < y < 1) presents a difficult and unresolved task. The systematic 
way of investigation of the intermediate cases (y ~ 1) is to use the Euler-Maclaurin 
summation formula and take into account the corrections in inverse powers of y(k) to the 
continuum limit. Such a treatment requires a numerical analysis. Alternatively, one may 
perform the RG studies by an integration in noninteger dimensionalities. In this Section 
we shall discuss the advantages and the disadvantages of this method. We shall show that 
it can be applied as an interpolation between the limiting cases j/<1 and |/> 1 only to 
specific theoretical schemes such as the RG. 
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4.1 5— integration 



Let us consider the integral 

y d-fc cthjy(j0/2] 
Mr,b)-L J (27r)d 2y(fc) , (12) 

which follows from (JHJ) and (JHJ. The integrand Si(fc, r) exhibits a single-power behaviour 
[~ y a (k)} only for <C 1 and ^> 1. The existence of a leading power dependence 
of the integrals A m+ \ = —(dA m /dr) on y = (Lq^) and the irrelevanve of the correction 
terms lead to a simple structure of the RG equations and, hence, to their scale invariant 
solutions, which are important for the description of the critical behaviour. The problem 
is to achieve such solutions for the intermediate cases of y ~ 1, too. 

It is impossible to construct an exact integral counterpart of the Eq. (J12)) with a power law 
behaviour with respect to y and for this reason, here we shall consider an approximate 
solution of the problem. We shall substitute the sum Si(k,r) with the 5-dimensional 
integral 

-5-1 f d s x 1 



S>1 ~ L ° J (2ny (x 2 + P + r) • (13) 
Accordingly, the integral will be substituted by the double (5, d) -dimensional integral 

\< i u\ t &— i f d k r d x 1 . . 

A i^b) = L J —J (2 ^ 5 ^ - - - r) . (14) 

The alternative is to substitute the integral f)12j) with the (d + <5)-dimensional integral 

^(r,6) = Lg-/ (2 ^ fa2 + r) , (15) 

where 5 G [0, 1]. Certainly, these substitutions are not exact counterparts of the original 
quantities and their utility in our attempts to present a reliable interpolation between the 
limiting cases should be justified. 

The integral (|T5j) can be deduced after the conjecture that the wave vector component 
k is a 8(< l)-dimensional (sub)vector, k = (k^; \i — 1, ...,5) and, accordingly, that the 
total wave vector q = {k^o'i K} and the volume Vd = L s Vd in the sum (|HJ) correspond to a 
(d + 5) -dimensional system [TT1 [TJl E] . The integrals A" and A" m+1 -^ = —(dA'^Jdr) have 
been used in preceding studies of quantum systems [12] and extended impurities [T3j for 
both 5 > 1 and 5 < 1. 

The integrals S[ and A', m+1 \ = —(dA' m /dr) given by (|TH|) and ()14|1 are defined with the 
help of another conjecture, namely, that one may perform a smooth interpolation between 
the integral values in the continuum limits (5 = 0) and 5 = 1 for the d- and D = (d+ 1)- 
dimensional cases, respectively, with the help of the formal rule 

^-/(S^f^" 1 ' (») 
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where K s = /it 5 l 2 T{8 /2). Using ©, © and (USD one immediately obtains the in- 
tegrals (pH?j) and dJ. The limit 5 — > in the last integral (in spherical coordinates) 
in (fTBj) should be taken with a special attention because of the divergency of the gamma 
function T(S/2). At first one should perform the integration over x of the integrand, 
say, x +0 f(x)/x, and then to take the limit 5^0. Usually, the integrands [~ f{x)] 
which appear by the perturbation series are such that no divergences arise in the final 
results for 5 ~ 0. This is confirmed by a direct calculation of the integrals ()13|). ()14|) and 
A' m+ i = -(dA'Jdr). 

For S — > and 5 — *■ 1, the integrals (fH?|) and (fT^j) exactly reproduce the results from (fTTj) 
and ()12j) for y — > and y — > oo, respectively. The same is valid for the integral (|15j) with 
respect to Eq. (|T^|) . So there are some grounds for the supposition that the intermediate 
states (y ~ 1) could be interpolated by the values < 5 < 1. 

4.2 Effective dimensionality of the fluctuation modes 

The coincidence of the results for the original integrals A m (r) in the limits y — > and 
y — > oo with the results from the integrals A' m and could be used as a basis of the 
supposition that there exists a continuous increasing function S(y) with the properties 
5(y — > 0) — > and 5(y — > oo) — > 1. This supposition presents the opportunity to 
intruduce a new dimensionality - the spatial dimensionality of the fluctuation modes 4>(q) 
given by 

D eS (y) = d + 5(y). (17) 

The Eq. (fTTj) is a straightforward generalization of the known (from previous FS stud- 
ies 0IH]) fact that the FS system abruptly changes its D-dimensional behaviour to the 
corresponding <i-dimensional behaviour when the thickness L is lowered to values less 
than £. According to the model (0), the system is represented by the field <j>{x). In this 
respect the notion for the effective length D e Q is not new but in this Section we consider 
it more explicitly within the generalized form given by Eq. (|17|) . Besides, in Section 5 
and 6 we shall introduce an e = (D e g — D^a) -expansion around the upper borderline 
effective dimensionality D^q. In general, this is a way to describe the RG scaling in 
terms of the fundamental ratio y = (L /£)- Unfortunately, the present investigation does 
not give an opportunity to obtain the function 8(y) and, hence, D c ^{y). The reason is in 
the pecuriality of the approach based on the 5-integrations introduced in Section 4.1. 

4.3 Validity 

Although the formal difference in their definitions, the integrals A' m (r) and A" m lead to 
the same results in many practical calculations. An example of such a calculation will be 
presented in this Section; for other examples, see Sections 5 and 6. The problem is that 
the values of the original sum (fTTj) and integral (fT2j) do not coincide with the values of the 
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corresponding integrals in noninteger dimensinalities (5^0). Therefore, the application 
of the integration in noninteger dimensionalities to concrete problems requires a special 
attention. The relaibility of such applications should be justified for any particular case. 
Despite the numerous RG studies which have been already performed with the help of 
the integrals A" m , the question about the limitations of the corresponding results has not 
been considered. In order to justify our RG investigations in Sections 5 and 6, we shall 
consider this problem. Besides, we shall demonstrate the degree of the approximation by 
compating the original sum S\(k,y) and integral Ai(y), denoting them as functions of y, 
with the corresponding integrals S[(k,y), A[(y), and A"(y). 

By calculating the integral (fTH|) for S[ and by comparing the result with Si from (JTTJ) we 
obtain 

th[y(k)/2]^g 5 [y(k)] (18) 

with 

1-5 

(19) 

where A(5) = n s / 2 /T(l — 5/2); hereafter the symbol "<->•" will denote a "comparison" 
and nothing else. The comparison (JTSJl can be done for any k < A. For k = 0, we have 
y(Q) = y, for k£ = 1, y(l/£) = v2, and for a third wave number k£ = v3, which also 
has an essential contribution in the integrals over k, we have y{V3/0 — 1- As we shall 
see our approximations are not precise enough in order to distinguish between the fits 
of Si and S[ for different wave numbers. Besides, there is a strong argument that the 
most important value of k at which we should make the comparison ()18|) is k = 0. This 
value corresponds to the uniform mode <p{k = 0) describing the spontaneous symmetry 
breaking in the c/-dimensional subsystem. 

By setting k = in (fTHj) and ()19|) we have 

th(y/2)^A(5)^ 5 . (20) 

The l.h.s and the r.h.s. of ()20|) are depicted in Fig. 1. From one side it is obvious that 
S[(0,y) is a good approximation to 5i(0, y) in quite broad intervals of values of y, for 
example: < y{k) < 1 and y{k) > 4. From the other side, it becomes evident that the 
values of 5 that give the best fit of ^'(O, y) to S(0, y) are: 5 = - for the case y < 1, and 
5=1- for the case y > 1. This means that, within the present consideration, the function 
5(y) can be approximated with zero for all y < 1 corresponding to a relatively good fit of 
the curves th(y/2) and gs(y) and that, owing to the same arguments, 5(y) ~ 1 for y ^> 2. 
In a broad region of values y around y = 2 the approximation of 5(0, y) ~ S'(0, y) breaks 
down and we could not deduce any reliable conclision about the exact critical value y c 
of y at which the FSC occurs. According to the present picture one may speculate that 
this value is probably y c = 2 whereas the intuitively appealing value is y = 1. The value 
y c = 2 comes from the factor (1/2) in the front of y{k) in Eqs. (JTSJ), (HHJ), and (|2T)j). The 
picture outlined from the comparison (|20j) is valid for the comparison denoted in (J19j) . 



g s [y(k)) = A(6) 



y(k) 
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In this case one should change y with y(k). For k = (v3/£), y c = 1, and the points of 
intersection of g^-lines will be located around the coordinate y c — 1. 

The calculculations within the present approach yield that the locations of the minimal 
values of the difference between the r.h.s. and the l.h.s. of (j2*U|) are 5 = for all y < 2 and 
5 = 1 for all y > 1, namely, that 5(y) coincides with the B-function Q(y — 2). Certainly 
this qualitatively wrong result is due to the obvious fact, see Fig. 1, that in a broad region 
around the point y c = 2, the approximation of 5*1(0, y) to S[(0, y) is not valid. It can be 
however concluded from such rough considerations that the real function 5(y) will have 
a steep increase from very small values 5 ~ up to 5 ~ 1 in a relatively close vicinity 
( 1 2/ — 2/c | ~ y c ) of the real critical point y c . 

The comparison of the integral Ai(r) with A[(r) and A"(r) is more difficult. In order to 
avoid the cutoff (A-) dependence of the results and simultaneously to avoid irrelevant to 
our problem ultraviolet divergences we shall set A = oo and consider the differences of 
type AAi(y) = [^4i(0) — Ai(y)] rather than the integrals themselves. Using Eq. (|T2^) we 
obtain 

AA^y) = -K d Ll d I d {y) (21) 



with 

where z - 
and flUJ): 

where 



h{y) 



dz z 



d-l 



cth(z/2) cth(v/PTF/2) 
z y/z 2 + y 2 



(22) 



L k. The difference AA' x {y) = [^(O) - A[(y)] is obtained from Eqs. © 
AA' 1 (y) = A(5)K d L 1 - d J' d (5 } y) } (23) 



dz z 



d-l 



1 



r 2-<5 



Finally, from Eqs. 



(z 2 + y 2 y- s / 2 

and (fTo^) we have that the difference [^"(O) — A"(y)] is given by 

„d+<5-3 



(24) 



POO 2 U 

AA'l(y) = K d+s Ll- d y 2 / dz 

Jo z z + y z 



(25) 



The integrals (J22j) . (|24p and (|25|) have obvious ultraviolet and infrared divergences at the 
corresponding lower and upper borderline dimensionalities. In order to avoid unnecessary 
complications in our calculation we shall consider the case of real films (d = 2). 

For (d = 2) the Eqs. (J23j) becomes 

A(5) 



&A'{(y) 



y_ 

2ttL 1 5 



(26) 



Because of the obvious infrared divergence at 5 
(dAA'l/dy) instead of the difference AA'{ itself: 

dAA'Ky) A(5)y 



0, we shall consider the derivative 



<5-l 



dy 



1t\Lc 



(27) 
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The ultraviolet divergences in the two parts of the integral J' 2 {5,y) are exactly compen- 
sated each other and this integral takes the simple form J' 2 (5,y) = y 5 /S. Thus we obtain 
that the derivative of AA[(y) is equal to the derivative (|27jl. The integral dl>i(y)/dy can 
be represented in the form 

Bhia) _ .. r » («Ht/2) + 1 \ (28) 



dy J y \ t 2 2tsh i (t/2) < 



which directly yields the result cth(y/2). For d = 2, the derivative of the difference 
AAi(y) becomes 

8My) - 1 -cth (I) . (29) 



dy AttL \2, 

The correspondence (|20|) straingtforwardly follows from the comparison of the Eqs. (}2Tj) 
and pnjl. 

Thus we have shown that the derivatives of the integrals and are quite different 
the derivatives of the original integrals A m , in particular for y ~ y c . Having in mind that 
the values of all integrals coincide at the limiting points 6(y = 0) = and S(y — oo) = 1, 
the same conclusion is true for the integrals themselves. 

The demonstrated deviation of the integrals in noninteger dimensionalities from the initial 
integrals A m does not mean that the RG analysis based on them is unreliable for < 8 < 1 . 
The argument here is that the RG predictions about the critical behaviour follow from 
the RG transformation which reflects the structure and the symmetry of the Hamiltonian 
rather than from the values of the perturbation integrals. The latters determine the 
location of the fixed points (FPs) of the RG equations and, therefore, might be of interest 
only in problems of special interest as to the question of whether a particular stable FP 
is accessible by the RG flows or not. For such special questions which often arise in 
investigations of complex systems the wrong determination of the FP coordinates, given 
by the integrals A' m or A'^, may produce wrong conclusions. Hopefully, our RG analysis 
in the next two Sections do not come upon such problems and we can arrive at reliable 
predictions about the critical behaviour of impure films irrespectively on the incorrest 
predictions about the location of the FPs. 

In Sections 5 and 6 we shall perform the RG investigation with the help of the integrals 
A' m . This more difficult variant is chosen to demonstrate that the results do not depend 
on the particular scheme of calculation. 



5 RG analysis: short— range impurity correlations 

The conventional RG treatment of the FSC has been presented in details in preceding 
works [7J |H] , where the perturbation integrals A m have been calculated by a direct sum- 
mations over k Q ; see, e.g. Eqs. and (JTJJ). We shall show that the application of 
the integration in noninteger dimensionalities yields more information about the critical 
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properties of films than the standard RG treatment of FS systems. For this aim we shall 
use the double (5, d)-integration. 

The RG recursion relations can be derived in two equivalent ways: by introducing an initial 
rescaling of the wave vector component k and without such an initial rescaling [ZHHj- We 
choose the latter variant in which the rescaling of the wave number k or, equivalently, of 
the thickness Lq ~ (1/ko), will appear as a result of the RG transformation. The initial 
rescaling transformations are k\ = bki and 

^ )a (k ,b- 1 k) = b 1 -^ 2 ( p' a (k ,k) , (30) 

where b > 1 is the rescaling number. In our investigation within the one-loop approxi- 
mation the Fisher exponent rj is equal to zero. The calculations are carried out by the 
direct way [""J E] of the averaging of the random functions. 

Using the standard way of calculations we have derived the RG recursion relations for 
short-range correlated point impurities (JH) in the form 

L' = &"% , (31) 
r' = b 2 {r + [4(n + 2)u - AJA'^r, b)} , (32) 



u 



b 3 - d {u - [4(n + 8)w 2 - 6wA]A' 2 (0, b)} , (33) 

A' = b 3 - d {A - [8(n + 2)uA - 4A 2 ]A' 2 (0, b)} , (34) 

where the integrals A' m (r,b) are the continuum limits of the sums (JBJ) with the upper 
cutoff A = 1 and an lower cutoff ft -1 ; hereafter we shall omit the (') of these integrals. 
The calculation of the integrals to first order in the small parameter r at the upper 
borderline dimensionality [D^g — (d + 5)u — 4] yields: 

Ai(r,6) = A x (0,6) - A 2 (0.b)r (35) 

with 

A JO h\ = T. s - x t(K\ I - 

2-5 



A l (0M = L 5 - 1 r(5)( 1 —^-) (36) 



and 

A 2 (0, b) = L s - l r(S)\nb , (37) 

where 

r(5) = (1 - 6 -)2A(5)K^ s ■ (38) 

Further, we substitute the results ([3*3)1 - (|3*7)l in the relations f)32|) - (f3*4*)l and perform the 
change of the parameters u and A with u = Lq~ 1 t(5)u and A = Lq _1 t(5)A. In the new 
notations, the Eqs. (|3*2*jl - (f3*4*)l take the form 

r' = b 2 {r + [4(n + 2)u - A] f — — — - rlnfo J , (39) 
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u' = b A - d ~ s {u - [4(n + 8)u 2 - 6uA]ln6} , (40) 

A' = b 4 - d - s {A - [8(n + 2)mA - 4A 2 ]ln6} . (41) 

The relation (J3"T|) describes the trivial scale invariance of the thickness L . This triviality 
is a result of the equivalence of the short-range random correlations along all D spatial 
directions (isotropy of the disorder). As we shall see in Section 6, the extended impurities 
break this isotropy and the RG relation for L becomes nontrivial. The structure of 
Eqs. (jnU) - (0H) is the same as that of the known RG equations for infinite systems [T5] . 
The difference is in the parameter 5^0. We shall see in the next Section that this simple 
structure of the RG equations reflects the above mentioned isotropy in the distribution 
of the disorder. Note, that the factor 1/(2 — 5) in Eq. (|39|) is relevant only for the precise 
determination of the FP value of the parameter r.The same iz valid for the 5-dependence 
of the integrals - (p?7|). In Section 4.3 we have already mentioned that the critical 
behaviour represented by the critical exponents does not depend on the concrete values 
of the FP coordinates in the parameter space (r,u,A) or, alternatively, in (r, u, A). 

The (d, <5)-integration is a relatively clumsy variant of calculation and cannot be easily 
applied in RG calculations in higher loop approximations. But the same results are 
obtained by the (d + <5)-integration (JT5j) which can be straightforwardly extended to 
considerations in higer orders of the loop expansion. Thus one can directly obtain the 
two-loop results for the same problem. The structure of the RG equations is again the 
same as that for the impurity problem in the corresponding infinite systems ^3]. The 
simple difference is that one should change d with (d + S), as it is for the Eqs. - (|41|). 

Thus we arrive at the main conclusion that irrespectively to the considered order of 
the loop expansion, the results for the impure films are obtained from those for infinite 
systems by the simple substitution of e = (4 — d) with e — [4 — (d + 5)] = (4 — -D e f). 
The same is valid for the FP coordinates. Alternatively, by using the e-expansion one 
obtains the results for both the D-dimensional hyperslab (5 > 0) and the corresponding 
infinite system (5 = 0). For the case of slabs we use the e = (d + 5)-expansion as shown 
in Eqs. (|39|) - (}4*Tj) and for the infinite system we take the limit 5 — > 0, which means 
to perform an expansion in e = e = (4 — d) = (5 — D). The reason for this (5 — D)- 
expansion for the infinite system is in the FSC (d — > D) denoting the equivalence between 
the critical behaviour in .D-impure films and that in d = (D — l)-dimensional impure 
infinite systems. Alternatively, in case of infinite systems, one may take the limit 6 — > 1, 
D b q = (d + 5) — > D and, therefore, the expansion parameter is e = 4 — D. 

Within the present RG analysis we have shown that the FSC can be described as a smooth 
change of the behaviour of the system with the variation of 8 from S = to 8 — > 1. Bearing 
in mind our discussion in Section 4, the variations of S between these limiting values mean 
a corresponding variation of the ratio y = (£o/0 from y <C 1 to y ^> 1. Although the 
smooth variation of 5 in the RG equations in the broad interval (0, 1), the FS crossover 
is expected to occur mainly in the interval \y — y c \ ~ y c , where the y c ~ 1 is the critical 
ratio. 
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In general, the critical behaviour of the hyperslab will correspond to the MF description 
for geometrical dimensionalities D > 5 — 5. As our RG results are valid above the lower 
effective borderline dimensionality (d + 8)l = 2, the nontrivial impure (for n < 4) or 
pure (for n > 4) stable critical behaviour will occur for geometrical dimensionalities 
(3 — 5) < D < (5 — 5). We shall not derive and discuss the FPs and the related critical 
and stability exponents of the impure films because we have demonstrated that these 
quantities and, therefore, the critical behaviour of the film as a whole, are obtained by 
setting e = (4 — d) — > e = (e — 8) in the known results for the respective quantities of 
impure infinite systems with the same type of disorder |13j . 

For real slabs (D = 3) we find that the nontrivial critical behaviour will appear for 

5 > 0. This implies L ^> a (quasi-2D films) and, moreover, the real films should be 
thick enough to ensure a parameter 5 > 0. In view of the discussion in Section 4, this 
probably corresponds to the ratio (Lq/£) of order of the unity. Denoting e = (e — 5), where 
e = (4 — d), we see that the results initially obtained for e< 1 can be extrapolated to real 
films (d = 2) provided we extend the values of our small parameter to 1 < e = (2 — 5) < 2. 
This is an advantage with respect to the simple extension of e = (4 — d) to e = 2 for 2D- 
films. There are no problems for the extrapolations of the e-results because the RG 
analysis does not demonstrate any peculiarities of the RG analysis at these relatively low 
dimensionalities such as: the appearance of new FPs, any runway of the knoun (at e <C 1) 
FPs, a qualitative change of the stability properties of the FPs or, a change of the location 
of the FPs in domains of the parameter space, where they could be unaccessible for the 
RG flows. 

From one side, we have shown, that the scope of the RG investigations with the help of 
field theretical models does not include 2D-dimensional films (mono-layers) or extremely 
thin films consisting of several atomic layers (L Q ~ a). From the other side, it has been 
completely clarified that the present approach can be undoubtedly applied to real films 
(D = 3). 

6 RG for extended impurities 
6.1 RG equations 

The RG equations for the case the extended impurities and inhomogeneities described by 
the Eqs. (J3J)— (|7J) are derived in the way outlined in the preceding Section. In order to 
simplify the notations we introduce the new variables v = Lq^ 1 t(5)u and /x = iQ-^A. 
Then, the RG relations can be written in the form 




(42) 
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b 2 < r + 4(n + 2)v 



- V 



1 - b- 2+5 



1-b- 2 
2-5 
b 8 - 1 N 
5 , 



— r In 6 



^ = 6 4 " d ^ - 8(n + 2)t^ In b + 4// 



'6*-l' 



(43) 

(44) 
(45) 



These RG relations have been obtained by the e = (4 — .D e g-)-expansion around the upper 
borderline dimensionality D^Z) = 4 (see Section 6). In contrast to the case in Section 6 
the scaling invariant solutions of these RG relations cannot be found for any 5 G [0, 1]. 
In order to obtain such solutions we must consider 5 <C 1 and to substitute the factor 
(b s -l)/5 with ln&. As our expansion parameter is e = [4 — (d + 5)] <C 1, this means 
that we must investigate the case when 5 <C and e = (4 — d) <S 1. Thus we have a single 
e-expansion with two small parameters: e = (4 — d) and 5. This variant of the theory is 
often reffered to as a "double (5, e)-expansion" [TQ1E1- The anisotropy of the random 
correlations "breaks" the single e-expansion to a double one. 

The //-contribution in the RG relation ()30|) for ko comes from a perturbation term shown 
diagrammatically in Fig. 2. In contast to the usual one-loop results, this self-energy 
contribution depends on the wave number k . The reason is that the broken line of the 
diagram in Fig. 3 does not carry external wave numbers k . This is a direct consequence 
of the fact that the random function ^p(k) does not depend on ko, i. e. of the infinite 
ranged correlations along the small site L . The /i-dependence in the RG equation for ko 
leads to a nontrivial scaling relation for the thickness L , namely, 



L' b 1+ ^ 2 



(46) 



which is quite similar to the scaling relations for the temperature T known from studies 
of quantum critical phenomena |17| . The scaling law (J34j) has the critical exponent 



1 + 



(47) 



where the asterisk (*) denotes any FP value of fi. This exponent is analogous to the 
dynamical critical exponent in the theory of dynamical critical phenomena in classical 
models with quenched disorder [TS] and quantum systems [TBI El ■ 

Using the relation ()46|) and 5 1, the Eqs. (J4*3*|) - (j45|) can be written in the form 



r' = b 2 1 r + [4(n + 2)v - //] 



V-2 



r In 6 



v 



I — h A-d-5(l+fi/2) 



\v - [4(n + 8)v 2 - 6vfi] lnfe} 



(4f 



(49) 
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fi' = b 4 - d - [8(n + 2)vfi - 4/x 2 ] In b} . (50) 

The extra-factor 6 -< W 2 j n (|49j) i s not essential for the RG analysis to this (e 1 -) order 
of the theory because fiS ~ e5 ~ e 2 . 

The RG relations (pSJl - (j5U|l yield new 5-corrections to the relevant physical quantities. 
They present the difference between the critical properties of the impure film and the 
corresponding infinite impure system. 



6.2 FPs and critical exponents 

The RG relations (f4*KJ) - (|5())l have four FPs: the Gaussian FP (hereafter referred to as 
GFP) with coordinates (ug = A*g = 0), the so-called "unphysical" FP (UFP), vy = 0, 
fi V = - e /4, the Heisenberg FP (HFP), 

6 3 

Vh = 77 — r~fiT ' Vh = 0, (51) 

4{n + 8) 

and the random FP (RFP), 

e + 25 (4 - n)e + 2(n + 2)5 

Vr = Tat TT ' ^ = a? n ' ( 52 ) 

16(n — 1) 8(n — 1) 

All these FPs are known by the work of T. C. Lubensky ^Hj- Here the FPs require a 
further investigation because of the FS effect (5 > 0). 

As usual, the GFP always exists and is stable for d > 4. The UFP has no physical 
meaning for e > 0, namely, for d < 4, because the parameter /i must always be nonnegative 
(fj, ~ A). For d > 4, the UFP is physical {p v > 0) but unstable. For 4 < d < (4 + 25), 
the UFP has an instability towards fi, whereas the same FP has a double instability 
(towards both v and [/,) for d > (4 + 25). Further, we shall concentrate our attention on 
the GFP, HFP and RFP. In order to analyze their properties we must obtain the critical 
and stability exponents. 

The stability exponents of the GFP are = e — 5 and Aj^ = e. These values show 
that the GFP is stable only for d > 4, where it describes an usual (free of disorder and 
fluctuation interactions) Gaussian behaviour. The stability exponents of the HFP are 

Af ) = 5 - e , (53) 

and 

Xf ) = < 4 - t + 2( " + 2M ■ (54) 
(n + 8) v ' 

( R) 

The stability exponents Ai 2 of the RFP are given by 

A S? = - a , 1 Tn fSne + 2(4 - n)«J =F V^l , (55) 
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where 



6 = (5n - 8) 2 e 2 - 4(15n 2 + 24n - 48)£ 2 

- 12(n 2 + 12n - 16)e<J . (56) 

The static critical exponents describing the critical behaviour of the system are given by 
rj = and the value of the correlation length exponent v. For the GFP we have vq = 1/2, 
for HFP, 

1 (n + 2) . , r _. 
^ = - + ^7 ^-(e-5), (57) 

2 4n + 8 V ; ' K J 



and for the RFP, 



1 3ne + 2(n + 2)5 



For 5 = and, hence, D e g = d, the Eqs. (I53 )l - (J55| ) yield the familiar results for d- 
dimensional infinite impure systems ^3] ■ For the GFP and HFP we have Z4 = 1, whereas 
for the UFP and RFP z ^ 1. For the RFP Egs. (g7j) and (J32D yield 



6.3 Stability properties of the FPs 

Consider the stabilyty properties of the HFP. The requirement of stability Aj^ < leads 
to the following inequality 

4(4 - d + 5) < (4 - d - 25)n . (60) 

When we solve this inequality with respect to n we should have in mind that the quantities 
(4 — d) = e ~ 5 are small and that is why it is not convenient to divide the ineq. (J6Uj) by 
the small factor (4 — d — 25). For this reason, let us consider first the case 5 = 0, where 
the inequality (|60|) is valid for n > 4, provided d < 4. Further, we consider n = 4 and 
for this value we easily find that the inequality (|6U|) cannot be satisfied for any 5 > 0. 
Therefore, for n = 4, the HFP is unstable towards the parameter fi, i. e. towards the 
disorder effects. For infinite systems (5 = 0) with n = 4, the HFP has a marginal stability 
(Aj^ = 0). The problem of whether the HFP is stable for n = 4 and 5 = has been 
widely discussed in preceding studies [12] in high orders in e. In particular, this topic has 
been extensively investigated by I. D. Lawrie et al (see Ref.[12]) for the case of impure 
systems with a cubic anisotropy. It has been shown that the higher orders in the loop 
expansion do not give the initially expected reliability in the treatment of the stability 
properties of the HFP near and at the value n = 4 of the symmetry index n. The reason 
is not in the specific features of the model but is a general disadvantage of the RG studies 
of systems with competing effects. Therefore, the same problem cannot be reliably solved 
in the present case, too. 
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The inequality ()6U|) is used together with A^ < and d > (2 — 5) in order to demonstrate 
that the HFP is unstable for all n < 4. For n > 4, the stability requirements A^ < 0, 
< and d > (2 — <5) yield the domain of stability 

2 - 5 < d < d H , (61) 

where 

<-» 

In fact, the conditions (|5T|) are satisfied for all 5 G [0, 1] provided n > 16, and for 
< 5 < do with 5q = 2{n — 4)/(n + 8) provided 4 < n < 16. Increasing the value of 5, the 
interval of variations of d decreases and for 5 = max(^, 1) the width of the <i-domain of 
stability becomes zero. 

The stability properties of the RFP can be investigated only for n ^ 1, because of the 
degeneration [T3| of the RG equations at n = 1. This degeneration leads to a special 
critical behaviour described by complex stability exponents ^3] and, hance, by oscillatory 
corrections to the main scaling laws ^3]- The investigation of the "oscillatory" critical 
behaviour is made in higher orders of e and 5. The values < n < 1 are of an academic 
interest only. We shall consider the case n > 1. 

Using the physical requirement fj,R > 0, one obtains an inequality inverse to (|50|l . Further, 
we must distinguish between two cases of stability of the RFP: (R) real negative exponents 
Xi2 < 0, and (C) complex exponets X[ 2 with negative real parts. For the case (R), using 
Eqs. flB3j) - (j5^j) . the requirement A^ < 0, d > (2 — 5) and the inequality inverse to (fBH|). 
we determine the following domain Rr of stability in the (d, n, 5) space: 

2 - 5 < d < (4 + 25) , (63) 
(4-n)(4-d) +2(n + 2)5 > , (64) 

3n(4 - d) + 2(4 - n)8 > , (65) 

and O > 0. For the case (ii) of complex stability exponents A^ 2 , the domain of stability 
Rc is defined by the inequalities (JMJ), (|63|) . 6 < 0, and d > (2 — 5). 

The criteria of stability can be investigated numerically for all allowed values of d, 5 and 
n^l. Here we shall restrict the analysis to those values of n which might be of importance 
for real systems or for the explanation of the properties of the model. The most important 
theoretical problem is the comparison of the domains of stabilities of the G, H and R FPS 
for values of n — 2, 3, 4, ... 

( R) 

For n = 4, one easily obtains the picture in Fig. 3. The RFP with real exponents X\ 2 
is stable in the domain R R defined by < 5 < <5 4 and (2 — d) < d < d R (4), where 
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<5 4 = 2/(1 + 2\/3) and c£r(4) = 4 — 2(1 + y/3). In Fig. 3, the point a has coordinates 
(S a = 5 4 ~ 0.45, d a = 1.55), the point 6 has coordinates (1, 1), and the point c is given 
by (1,4). The domain Rc in Fig. 3 denotes the stability region of the random critical 
behaviour with complex stability exponents. This domain is given by the inequalities 
d R (A) < d < 4 provided 5 < 5 4 , and by (2 - 5) < d < 4 for 5 4 < 5 < 1. 

For 1 < n < 4 the stability criteria for R in the case of real stability exponents are > 
and 

2 - 5 < d < d R , (66) 

where 

The domain i?c of complex stability exponents is given by the inequalities < and (|66p. 

The domains of stability .Rr and i?c for 1 < n < 4 are quite similar to those in the case 
n — 4 depicted in Fig. 3, but the coordinates of the points a and c in the (5, d) plane vary 
with n. The coordinates (5 a , d a ) and (d c , 5 C = 1) of the points a and c are given in Table 
1 for several values of n. The value 5 a decreases from 0.45 to zero with the decrease of n 
from n ~ 4 to n ~ 8/5. For n < 8/5, S a increases with the decrease of n. 



TABLE 1. Values of 5 a , d a , and d c 



n 


7/5 


8/5 


9/5 


2 


3 


4 




1.06 




0.02 


0.05 


0.26 


0.45 


d a 


1.94 




1.98 


1.95 


1.74 


1.55 


d c 


5.24 


5 


4.81 


4.67 


4.22 


4 



The value n = 8/5 is quite special because for this value the function 0(n, e,5) given 
by Eq. (|56|) depends linearly on e. Within our one-loop approximation, the domain Rr 
at n = 8/5 does not exist at all and the domain Rc is expanded up to < 5 < 1 and 
2 — 5 < d < 4 + 5. In our case, the effective values of 5 a (8/5) and d a (S/5) are and 2, 
respectively. It has been shown by I. D. Lawrie et al ^H] that the same value n = 8/5 
is the reason for the peculiar critical behaviour of infinite impure systems (5 = 0) with 
a cubic anisotropy. So, one may expect that in our case a similar peculiar behaviour at 
n = 8/5 with singularities of the e-expansion and complex exponents will occur. 

It is seen from Table 1 that the coordinate d c (n) in Fig. 3 is larger than d c (4) = 4 for all 
< n < 4. This means that for all 1 < n < 4, there is a triangle domain for certain values 
d > 4, where both the RFP and the GFP are stable. In this region there is a competition 
between the pure Gaussian behaviour represented by the GFP and the impure behaviour 
represented by RFP. The outcome of this competition should depend on the strength of 
the disorder effect. 

For n > 4, both the RFP and the HFP are stable for certain values of d and S. The 
stability domains Rr and Rc can be determined with the help of the inequalities (jo^|) - 
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(JBSJ). For 5 < 5q, the domain R R is confined by the inequalities > 0, and 

d H < d <d R (68) 
whereas for 5 > So, the same domain is described by the inequalities 6 > and 

2 - 5 < d < d R (69) 

The domain R c is defined by the same inequalities (jSHj) and (jSHj) but for 6 < 0. 

A typical picture of the stability domains in this case is given in Fig. 4 for n = 6, 
where the points h and a have coordinates: 5h = S (6) = 0.29, dh = 1.61, 5 a = 0.74 
and d a = 1.26. The HFP is stable for relatively small values of 5, whereas the stability 
of the RFP dominates for large values of 5. There is a domain (4 — 0.225) < d < 4 of 
dimensionalities d ~ 4 for which both the RFP and the HFP are unstable. This domain is 
defined by the location of the point c with coordinates S c = 1, d c (6) = 3.78; see the shaded 
region in Fig. 4. Such domains of instability exist for all symmetry indices n > 4 and 
5 > 0. They can be described by the inequality d'(n,5) < d < 4 which can be obtained 
by the RG analysis. For n = 6 and 5=1, d'(n, 5) = d c (6) as shown in Fig. 4. For <5 — > 0, 
d r (n, 5) — > 4 and for 6 — > 1, d'(n, 6) is lowered up to a minimal value d'(oo, 1) = 3.33 (See 
also Fig. 5). 

As the GFP is stable only for d > 4, the shaded domain in Fig. 5 remains the region of 
a total instability of the system. This total lack of stable FPs for all n > 4, 5 > and 
certain d ~ 4 can be interpreted as an indication for a fluctuation-driven impure tricritical 
phenomenon followed by a first-order phase transition. The FS system (5 > 0) is unstable 
towards the disorder (the parameter /i) and, hence, the reason for the appearance of a 
fluctuation-driven phase transition of first-order is in the simultaneous effect of the FS 
and the extended impurities. 

In order to justify the prediction for the existence of a tricritical point and a fluctuation (or 
disorder)-driven phase transition of first order, we remind that the fluctuation interactions 
represented by the parameter u are relatively small near four dimensionalities (d ~ 4). 
Under such circumstances, the disorder effects may alter the sign of the effectivelly small 
parameter u from u > (a second-order transition) to u < (a first-order transition). 
At u = 0, a tricritical point should appear. Note, that these effects are impossible for 
infinite impure systems and in pure FS systems. The change of the order of the phase 
transition in the present case is a result of the competition between fluctuation, disorder 
and FS effects. The effective fluctuation model for the description of this change of the 
order of the phase transition can be constructed by adding an additional Uq4> 6 term to 
the Hamiltonian (Q). 

The stability domain of the HFP vanishes when the symmetry index n decreases from 6 
to 4. Increasing the number n to n = 16, the stability region of the HFP expands up to 
the triangle 264 (see the dotted line 46 in Fig. 4). The further increase of the values of n 
above n = 16 leads to an enlargement of the stability region of the HFP so that the points 
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h and a already lay on the vertical line 5 = 1. For n — > oo (the limit of the spherical 
model [Hj) the point h has coordinates (1,2) on the vertical line 5 = 1, as shown in Fig. 
5. The Rr is extremely small but persists (see the narrow wing Ma with d a = 2.19). The 
domain Rq is shown in Fig. 5 by the triangle a4c, where d c = 3.33. In Figs. 4-6, the 
point b has the same coordinates (1,1). 

7 Summary and validity of the results 

We have obtained and analyzed in details the RG recursion relations for FS systems of 
slab geometry described by the n-component LG model containing quenched impurities. 
Two models of quenched impurities have been considered. The results demonstrate the 
outcome of the competition between the effects of fluctuation interactions, the finite size, 
and the disorder. 

We have used an e = (4 — .D c g-)-expansion defined by the effective spatial dimensionality 
S = (D e g — d) of the fluctuation modes of the order parameter. Our approach is based 
on an approximate substitution of the lattice summation with an integration in dimen- 
sionalities less than unity. The limitations of this approach have been considered. The 
approach is more general that the familiar RG considerations based on the e = (4 — d)- 
expansion. The fundamental relationship between the effective dimensionality _D e g- and 
the ratio y = (Lq/£) of the characteristic lengths of the film has been introduced and 
discussed. 

The most interesting cases of critical phenomena are those for 5^1 (thin films) and 
5 ^> 1 (thick films, when the film behaves as an almost-infinite system). It has been 
demonstrated that our approach based on an integration in dimensionalities less than 
unity is reliable for the description of these cases, in particular, for real films (D = 3). 

The method introduced in this paper has been used to demostrate the validity of the FSC 
in films with homogeneously distributed quenched impurities with short-range random 
correlations. It has been shown that the proof of the FSC in these systems can be easily 
proven and straightforwardly expanded to any order in the loop expansion. The reason 
is that the impurity correlations are spatially isotropic and the sum (e + 5) can be taken 
as the single expansion parameter e. 

The extended impurities break the spatial isotropy of the disordered system, and the 
RG investigation yields a dependence of the critical exponents on two small parameters: 
e = (4 — d) = (5 — D) and 5. This leads to another picture of the critical behaviour which 
will be summarized in the remainder of this Section. 

The dependence of the FP coordinates and the critical exponents on the small parameters 
(e, 5) and the symmetry index n has been calculated and analyzed with respect to the 
stability properties of the FPs. In this way we have established the possible types of critical 
behaviour for different values of 5, the geometric dimensionality D of the hyperslab, and 
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n. In order to describe the random critical behaviour in films with extended impuirities 
we have introduce and calculated a new critical exponent z&. The analogy between this 
exponent and the dynamical exponent in quantum systems has been emphasized. The 
results are valid for effective dimensionalities D e g > = 2 - the lower borderline 

dimensionality. 

For real films (D = 3) with extended impurities our analysis is valid for the so-called 
quasi-2.D (thin films), where 5 > 0. So, the results from our RG analysis performed for 
small values of e and 5 can be exptrapolated for predictions of the critical behaviour of 
real impure films (D — 3, d — 2) with a finine thickness Lq ^> a, where the S(y) > 0. 
We have shown that, in the strict mathematical sense, the latter condition is always valid 
for y > 0, i. e. when the system is outside the critical point. But in real expriments 
the conditions for the validity of our consideration are somewhat different. The function 
5(y) ~ up to values of y ~ 1 and, therefore, in experimental conditions (£ ^> a) the 
thickness L ~ £ should be much larger than the interparticle distance a. As far as 
the behaviour in the critical region (£ ^> a) is concerned the condition y > will be 
easily satisfied for films with L ^> a. This is the case of interest for experiments. The 
case of exactly 2D films (Lq = a) corresponds to the lower borderline dimensionality 
and, therefore, is beyond the scope of our investigation. The invalidity of the results at 
the lower borderline dimensionality is a generic disadvantage of the field-theoretical RG 
methods rather than a result of our particular approach to the problem. The 2D case can 
be effectively achieved for any film in the asymptotic vicinity of the critical point, where 
£ — > oo and, hence, the ratio y — > for any finite thickness Lq. However, this situation is 
not of interest for the most part of real experiments. 

While the condition D c g > D^ = 2 guarantees the reliability of the qualitative pre- 
dictions of our RG analysis, the extrapolation of the results obtained near the upper 
borderline dimensionality D^q = (4 — d — 5) to the real dimensionality (d = 2) requires 
an extension of the values of the small parameter e to values e ~ (2 — 5) < 2. This is 
the usual way of extrapolation of the RG results to real systems. It is believed the small 
parameter e should be set equal to (2 — 5) in order to predict the values of the critical 
exponents for the real case (8 <^ 1) and that in this case, the one-loop results do not 
give quantitatively correct predictions. Although this is generally true, the qualitative 
predictions from the one-loop results are reliable, in particular, for the solution of im- 
portant problems like the possible types of critical behaviour and the conditions under 
which a critical behaviour can occur. In this article we have been mainly involved in these 
type of problems and the results are well grounded within the one-loop approximation. 
The practice of the numerous applications of the RG to complex systems shows that the 
one-loop results about the general picture of the critical phenomena remain valid within 
the framework of consideration by higher loop approximations. For real films e = (2 — 5) 
should take values e ~ 2 for 5 1 and e ~ 1 for 5 ~ 1 (thick films). Therefore, the analy- 
sis carried out in Section 6 is reliable for the case of real films and their critical behaviour 
is given by Figs. 3 - 5 on the line d — 2. In addition to this conventional point of view 
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we shall stress that the results for small values of e = (4 — d) can also be extrapolated to 
real films but for those parts of the critical region which are not asymptotically close to 
the critical point. Thus the fluctuation-driven change of the order of the phase transition 
near the upper borderline dimensionality (4 + 5) can be extrapolated to real effective 
dimensionality (d + 5) , provided one is interested in fluctuation phenomena which occur 
outside the asymptotic vicinity of the critical point. 

A brief summary of other results for the critical behaviour in real dimensionalities (d = 
2,0 < 5 < 1) can be outlined with the help of Figs. 3-5. This critical behaviour is 
described by the HFP and the RFP. For 1 < n < 4, only the random critical behaviour is 
present. For relatively thick films this critical behaviour will exhibit oscilatory corrections 
to the scaling laws described by complex stability exponents, whereas the ramiliar random 
critical behaviour will occur in relatively thin films. For n > 4, the pure critical behaviour 
represented by the properties of the HFP will occur for relativelu thin films whereas the 
random critical behaviour in its two variants of real and complex stability exponents will 
occur in thick films. In the spherical limit (n = oo) the pure critical behaviour is stable 
for all real films and the disorder effects are irrelevant for the critical properties. The 
singular behaviour of the films for the symmetry indices n = 1 and n = 8/5 has been 
explained in Section 6. 

Finally, we note, that the present investigation reveals essential differences in the critical 
behaviour of infinite and FS systems with extended impurities. The finite thickness L 
of the film leads to: an essential dependence of the critical exponents on the noninteger 
dimensional portion 5, a fluctuation-driven change of the order of the phase transition for 
all n > 4 and d ~ 4 (the almost-MF region), the competition between the pure (Gaussian) 
and impure critical regimes for all 1 < n < 4 and 2 < d < 4, and the dominating role 
of the impure (random) critical behaviour with oscillatory corrections to the scaling laws 
for all relatively thick films. 

Acknowledgments: 

L. Craco and D. I. Uzunov gratefully acknowledge the hospitality of the MPI-PKS (Dres- 
den), Aussenstelle Stuttgart and, in particular, the kind support by Prof. Peter Fulde. 
L. Craco thanks the support by the Brazilian Agency Conselho Nacional de Desenvolvi- 
mento Cientifico e Tecnologico (CNPq). A research grant Ph560 of NFS (Sofia) is also 
acknowledged. 



24 



References 



[1] See, e. g., K. Binder, in: Phase Transitions and Critical Phenomena, Vol. 8, ed. by C. 
Domb and J. L. Lebowitz (Academic Press, London, 1983) p. 1; H. W. Diehl, ibid, 
Vol. 10, p. 76. 

[2] J. J. Pearson, Phys. Rev. 138 (1965) A213; D. J. Mills, Phys. Rev. B 3 (1971) 3887; 
R. A. Weiner, Phys. Rev. B 8, (1973) 4427; R. Pandit and M. Wortis, Phys. Rev. B 
25 (1982) 3226. 

[3] D. I. Uzunov and M. Suzuki, Physica A 204 (1994) 702. 

[4] M. I. Kaganov and A. N. Omel'yanchuk, Zh. Eksp. Teor. Fiz. 61 (1971) 1679 [Sov. 
Phys. - JETP 34 (1972) 895]; T. C. Lubensky and M. H. Rubin, Phys. Rev. B 12 
(1975) 3885. 

[5] M. E. Fisher, in Critical Phenomena, Proceedings of the 51th Enrico Fermi Summer 
School, Varenna, Italy; ed. by M. S. Green (Academic Press, New York, 1971) p. 1; 
M. E. Fisher and M. N. Barder, Phys. Rev. Lett 28 (1972) 1512; 

[6] For a review, see, e.g., M. N. Barber, in: Phase Transitions and Critical Phenomena, 
Vol. 8, ed. by C. Domb and J. L. Lebowitz (Academic Press, London, 1983) p. 145. 

[7] See, e.g., J. A. Bray and M. A. Moore, J. Phys. A 11 (1978) 715; I. D. Lawrie, J. 
Phys. C 11 (1978) 3857; A. M. Nemirovsky and K. F. Freed, J. Phys. A 19 (1986) 
591; H. Guo and D. Jasnow, Phys. Rev. B 35 (1987) 1846; D. O'Conner and C. R. 
Stephens, Nucl. Phys. B 360 (1991) 297. 

[8] D. I. Uzunov and M. Suzuki, Phisica A 216 (1995) 489. 

[9] D. I. Uzunov, Introduction to the Theory of Critical Phenomena (World Scientific, 
Singapore, 1993). 

[10] See, e.g., G. Grinstein, in Fundamental Problems in Statistical Mechanics VI, ed. by 
E. G. D. Cohen (North-Holland, Amsterdam, 1985) p. 147; M. McCoy and T. T. 
Wu, Phys. Rev. Lett. 21 (1968) 549; Phys. Rev. 176 (1968) 631; 188 (1969) 982; 188 
(1969) 1014. 

[11] L. De Cesare, L. Craco, I. Rabuffo, and D. I. Uzunov, J. Phys. Studies (1998), in 
press. 

[12] A. M. J. Schakel, Phys. Lett. A 224 (1997) 287; T. R. Kirkpatrick and D. Belitz, 
Phis Rev. Lett 76 (1996) 2571. 

[13] T. C. Lubensky, Phys. Rev. B 11 (1975) 3573; D. E. Khmelnitskii, Phys. Lett 67A 
(1978) 59; I. D. Lawrie, Y. T. Millev and D. I. Uzunov, J. Phys. A 20 (1987) 1599; 
See also S. Ma, Modern Theory of Critical Phenomena (Benjamin, London 1976). 



25 



[14] S. N. Dorogovtsev, Phys. Lett. 76A (1980) 169; D. Boyanovsky and J. L. Cardy, 
Phys. Rev. B 26 (1982) 154; 27 (1983) 6971(E); I. D. Lawrie and V. V. Prudnikov, 
J. Phys. C 17 (1983) 1655; A. Weinrib and B. I. Halperin, Phys. Rev. B 27 (1983) 
413; PB. Weichman and K. Kim, Phys. Rev. B 40 (1989) 813. 

[15] G. Grinstein, S. Ma and G. Mazenko, Phys. Rev. B 15 (1977) 258. 

[16] E. R. Korutcheva and D. I. Uzunov, Phys. Lett 106A (1984) 175; G. Busiello, L. Dc 
Cesare, and I. Rabutto, Phys. Lett. 102A (1984) 41. 

[17] J. A. Hertz, Phys. Rev. B 14 (1976) 1165; G. Busiello and L. De Cesare, Phys. Lett. 
77A (1980) 177; D. I. Uzunov, Phys. Lett. 87A (1981) 11; A. Caramico D'Auria, L. 
De Cesare, and I. Rabutto, Physica A 225 (1996) 363. 



26 



Figure captions 

Fig.l. First-order diagram for k' . The discontinuous line reprsents the disorder average 
[|(/?(A;)| 2 ]# = A and the continuous line represents the bare correlation function 

— * — * 

G (k , k + ki). 

Fig. 2. Stability domains of R for n = 4: R R (triangle 2a4) and R c (Aabc). 

Fig. 3. Stability domains of H (triangle 2M) and R for n > 4. Rr is given by the triangle 
Ma and Rc is the rectangle a4cb. The dotted line 64 marks the extension of the stability 
domain of H for n = 16. The shaded region is explained in the text. 

Fig. 4. Stability domains of H and R for n = oo (Hartree limit [10]). Rr and Rc are 
given by the triangles Ma and a4c, respectively. 
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